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ABSTRACT 


Avail  uiuv  Codr 


One-dimensional,  steady  flame  propagation  for  a  sequential,  two-step 

Xu 


X 


reaction  of  the  form  A  -*•  B  +  C  is  considered.  An  earlier  investigation  of 
the  problem  by  Kapila  and  Ludford  (Combustion  and  Flame  29 ,  p.  167  (1977)) 


determines  that  two  separated  flames  generally  exist  and  that  their  ordering 
is  fixed  by  the  ordering  of  the  (disparate)  magnitudes  of  the  activation 
energies.  The  present  work  shows  to  the  contrary  that  reversals  of  the  flame 
ordering  are  quite  possible,  but  that  this  is  a  subtle  effect  requiring 
attention  to  issues  which  are  usually  ignored  in  the  theory  of  single  flames. 


AMS  (MOS)  Subject  Classifications:  80A25,  34E05 

Key  Words:  two-step  sequential  reaction,  large  activation  energies, 
matched  asymptotic  expansion 
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SIGNIFICANCE  AND  EXPLANATION 


Asymptotic  methods  based  on  large  activation  energies  have  been 
established  over  the  last  decade  as  an  effective  technique  for  analyzing 
chemical  processes  in  combusion  which  are  described  by  Arrhenius  kinetics. 
However/  efforts  have  been  largely  restricted  to  simple  single-step  reactions, 
which  are  of  limited  practical  interest.  The  need  for  a  more  thorough 
treatment  of  multiple-step  reactions  has  been  recognized  for  some  time  and 
motivates  the  present  investigation  as  a  step  towards  narrowing  the  gap 
between  the  basic  theory  and  the  practical  needs  of  engineering  applications. 

The  contribution  of  this  report  to  the  analysis  of  multiple-step 
combustion  is  two-fold .  It  provides  a  careful  treatment  of  the  title  problem 
and,  more  importantly,  it  demonstrates  the  sometimes  subtle  complexities 
introduced  by  the  presence  of  more  than  one  reaction.  The  need  for  attention 
to  detail  is  exemplified  as  subtle  differences  in  analysis  are  shown  to  lead 
to  remarkably  different  conclusions. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


TWO-STEP  SEQUENTIAL  REACTIONS  FOR 
LARGE  ACTIVATION  ENERGIES-REVISITED 

H.  V.  McConnaughey  and  G.  S.  S.  Ludford 

INTRODUCTION 

This  discussion  is  concerned  with  steady,  one-dimensional  propagation  of 
the  two  laminar  flames  associated  with  the  two-step  sequential  reaction 
A  +  B  ♦  C.  Of  particular  interest  are  the  admissible  flame  configurations 
and  associated  burning  rates. 

Kapila  and  Ludford  [1]  investigated  this  problem  using  an  asymptotic 
analysis  based  on  the  limit  of  large  activation  energies,  E1  and  Ej.  It 
was  found  that  two  separated  flames  generally  exist,  but  that  the  flames  can 
also  merge.  The  order  in  which  the  separated  flames  occur  was  shown  to  be 
fixed  by  the  relative  sizes  of  the  activation  energies.  In  all  cases,  the 
burning  rate  was  determined  explicitly. 

This  report  presents  a  generalized  version  of  [1]  which  also  considers 
the  limits  E^  +  •  and  E2  +  «.  The  generalization  is  subtle,  however, 
involving  only  terms  which  are  OtE^)  or  OCE^),  and  is  likely  to  appear 
superfluous  at  first  glance.  Indeed,  it  reproduces  the  results  found 
previously  for  the  case  of  separated  flames  and  yields  burning-rate  formulas 
for  merged  flames  which  are  nearly  identical  to  those  in  [1].  Nevertheless, 
the  restriction  on  the  separated- flame  ordering  found  by  Kapila  and  Ludford  is 
eliminated. 

The  nature  and  significance  of  this  generalization  are  best  seen  in  the 
context  of  the  asymptotics,  which  are  therefore  included  in  this  discussion  to 
the  extent  necessary.  Treatment  of  separated  flames  is  omitted  since  it  does 


This  research  was  sponsored  in  part  by  the  United  States  Army  under  Contract 
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not  appreciably  modify  the  conclusions  derived  by  Kapila  and  Ludford  from 
their  investigation  of  separated  flames  (see  [2]).  An  analysis  of  merged 
flames  is  needed,  however,  and  hence  is  presented  here  along  with  pertinent 
results  for  separated  flames.  The  approach  in  [1]  is  followed,  thus  the 
reader  is  referred  to  that  work  for  any  necessary  clarification. 

It  should  be  emphasized  that  most  of  the  present  analysis  follows  that  of 
Kapila  and  Ludford,  although  extra  terms  may  appear  and  the  notation  may 
differ  slightly.  The  significant  deviations  from  their  work  are  explicitly 
stated . 

The  Mathematical  Problem 

Steady  plane  flame  propagation  for  a  sequential  two-step  reaction 
A  ♦  B  C  is  considered.  The  governing  dimensionless  equations  are 


where 


L(Y)  *=  , 

L(Z)  -  -W^  +  , 

L(T)  «  -QtW'  -  Q2V’2  , 


x  6  (-«,<»)  , 
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and  where 

Y  »  mass  fraction  of  species  A, 
Z  •  mass  fraction  of  species  B, 


T  -  dimensionless  temperature 


subscript  1  ~  pertaining  to  reaction  A  ♦  B, 
subscript  2  ~  pertaining  to  reaction  B  C, 

expfE^/T^)  *  =  Damkohler  number, 

ni 

*  E^  ,  where  nt  *  unspecified  constant, 

E^  «  dimensionless  activation  energy, 

T^  »  parameter  characterizing  magnitude  of  D^, 

M  «  burning  rate, 

Qa  -  dimensionless  heat  release,  Qi  +  C*2  ” 

(The  reader  is  referred  to  [1]  for  additional  details.)  The  boundary 
conditions  are 

X  ♦  -*s  Y  ♦  Y^,  Z  +  Z_^,  T  *  T^  |  (4) 

x  ♦  -H»:  Y  ♦  0,  Z  ♦  0,  T  ♦  T  .  (5) 
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Note  that  this  system  is  invariant  under  translations,  therefore  the  origin 
may  be  fixed  as  desired.  Also,  one  of  equations  (1)-(3)  may  be  eliminated  by 
integrating  the  combination  (1)  +  Q2{2)  +  (3)  subject  to  (4)  and  (5).  The 
resulting  identity: 

Y(x)  +  Q2Z(x)  +  T(x)  =  Y^  +  -  T.  (6) 

holds . 

The  objective  of  Kapila  and  Ludford's  investigation  is  to  determine  from 
( 1 )  — ( 6 )  the  variation  of  M  with  D1  and  D2,  the  other  parameters  being 
fixed.  This  is  accomplished  by  an  asymptotic  analysis  based  on  large 
activation  energies.  In  the  limits  E^  ♦  «  and  ♦  •,  the  chemical 
activity,  represented  by  the  nonlinear  terms  in  (1)-(3),  is  localized  in  thin 
zones  (flames)  where  the  temperature  is  close  to  T^  or  T2.  Outside  of 
these  zones,  the  linear  reactionless  form  of  equations  (1)-(3)  holds.  The 
resulting  "outer"  solutions  are  matched  at  the  flames  to  the  solutions  of  the 
(nonlinear)  equations  valid  inside  the  reaction  zones. 


Two  types  of  flames  are  possible:  1)  separated  flames,  corresponding  to 
cases  where  -  Tj  ■  0(1),  by  which  the  combustion  field  is  divided  into 
three  reactionless  regions,  and  2)  merged  flames,  which  separate  two 
chemically  inert  regions  and  occur  when  ■  OfE^1),  Ej  «  minCE^Ej). 

<v  *** 

For  convenience,  we  introduce  the  constants  T^  and  T^,  where 

T^  *  lim  T^  . 

E  .E  x® 

1  '2 

In  this  notation,  separated  flames  exist  when  T^  *  T^;  merged  flames  occur 
if  T,  «  T2. 

Kapila  and  Ludford  do  not  distinguish  between  T^  and  when  k  is 

the  index  associated  with  the  reaction  which  exhausts  all  reactants  (i.e. 
beyond  which  Y  -  2  *  0).  For  that  index,  T^  «  T—  and  in  [1],  it  is 
effectively  assumed  that  is  strictly  equal  to  Tw .  it  is  this  assumption 

which  leads  to  the  aforementioned  restriction  on  the  ordering  of  separated 
flames,  as  will  be  shown. 

Results  for  Separated  Flames  (T1  *  Tj) 

The  reaction  occurring  at  the  lower  temperature  precedes  (in  location) 

<v 

the  other  reaction.  For  T^  <  ,  the  flames  then  appear  as  shown  in  Figure 

/■v 

1  and  T2  equals  T^.  ihe  burning  rate  for  this  case  obtained  in  [1] 
satisfies 
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For  T  <  T  ,  the  flames  in  Figure  1  are  reversed,  T.  *  T  ,  and  M  is 
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Merged  Flames  (T1  -  T2> 

When  the  flames  are  merged,  the  combustion  field  is  divided  into  only  two 
chemically  inert  regions,  as  represented  in  Figure  2a.  The  origin  has  been 
located  within  the  flame  zone.  The  solution  of  the  reactionless  equations  for 


x  *  0  subject  to  (4)  and  (5)  is 

x  <  0;  Y  =  ( 1  -  eX)Y  +  veX,  Z  *  (1  -  eX)Z 

—00  —00 


ve 


Y  -  P2Z  , 


x  >  0:  Y  *  0,  Z  -  0,  T 


-1 


where  v  is  an  unspecified  quantity  of  order  E  Continuity  of  Y  and 
Z  to  leading  order  across  x  =  0  has  been  imposed;  the  origin  has  been  fixed 
so  that  Y  +  Z  is  independent  of  E^  and  Ej.  The  leading-order  flame 

temperature  T,  =  T„  must  equal  T  . 

12  00 

The  validity  of  the  above  outer  solution  is  contingent  upon  the  existence 
of  a  flame-zone  solution  or  "structure"  in  the  neighborhood  of  x  **  0  which 
effects  the  change  in  slope  there.  The  structure  analysis  below  is  seen  to 
fix  the  burning  rate  M. 

Equations  ( 1 )  and  ( 2 )  indicate  that  the  A  B  reaction  and  the  B  +  C 
reaction  become  active  when  T  -  T^  *  (ME^),  i  =  1,2.  The  small  parameters 

6.  and  $  which  gauge  the  thickness  of  the  flames  are  thus  chosen  so  that 
1  2 

6^  *  OfE^  ).  Although  the  two  reactions  occur  at  the  same  0(1)  location 
(x  *  0),  they  must  be  distinct  in  order  to  make  the  structure  problem 
analytically  tractable.  It  is  therefore  assumed  that  the  reaction  zone 
thicknesses  6^  and  6^  are  of  disparate  orders,  or  equivalently,  the 
activation  energies  are  of  disparate  orders. 

Consider  the  case  where  >>  E2,  so  that  <<  ^2’  7116  associate<J 

merged-flarae  configuration  is  illustrated  in  Figure  2b.  On  the  scale  of  the 
(broader)  B  ♦  C  reaction,  the  A  ♦  B  reaction  is  located  at  x/6„  *  p. 


where  it  appears  as  a  discontinuity.  All  chemical  activity  of  the  A  ♦  B 


L ( Y) s  0  ,  L(Z)=  0 


Y  =  0  ,  Z  =  0 


»  x 


x  s  0 


Figure  2a.  Appearance  of  merged  flames  on  the  0(1)  scale:  A 
single  discontinuity  at  x  =  0. 
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Figure  2b.  Merged-f lame  configuration  for  >>  E2:  A  -*■  B 
reaction  zone  appears  as  a  discontinuity  at  p0 
on  the  scale  of  the  broader  B  -*■  C  reaction. 


reaction  is  localized  and  completed  there;  the  reaction  is  frozen  for 
x/&2  *  P  <  P0  and  is  in  equilibrium  (reactant  A  is  exhausted)  for 
p  >  Pq.  The  B  C  reaction,  on  the  other  hand,  is  active  for  p  e 
It  is  only  partially  completed  for  p  <  p^  and  is  not  significantly  affected 
at  P0  by  the  A  ♦  B  reaction  there  since  the  latter  occurs  in  such  a  thin 
zone.  The  amount  of  reactant  B  remaining  at  p  =  p^,  and  all  B  produced 
at  pQ,  are  completely  converted  into  C  over  the  interval  (p0,»)  of  p. 
Thick -Flame  Structure 

In  the  B  C  reaction  zone,  the  variables  are  written: 
x  =  62p,  p  e  , 

T  =  T^  +  S2t(p)  +  o(fi2>  , 

Z  =  62z(P)  +  o(S2)  , 

f  S2y(p)  +  o(S2)  for  P  <  PQ 


Y  =  i 


(.  0  for  p  >  p0  , 

where  62  =  T«/E2'  =  "WP)  ~  C?2z(p)  for  p  <  pQ  and  t(p)  *  -Q2 z(p ) 

for  p  >  pQ.  The  parameters  T1  and  T2  may  similarly  be  expressed  as 

T1  *»  T^  +  ^2tl  +  °^2^'  T2  ~  T«»  +  ®2fc2  +  *>(62)  • 


(Kapila  and  Ludford  assume  t2  -  0.)  To  leading  order,  equations  (1)  and  (2) 
are  reduced  to 


P  <  PQ  8  “2  "  52D1M  2y(P  )exP^t{p  >  “  t1]E1/B2>  ,  (9) 

dp 

2  2 

+  «2D^M~2z(p)exp(t(p)  -  t2)  , 

dp  dp 

and 

2 

p  >  pQ  :  Y  *  0,  -  62D£M~2z(p)exp(t(p)  -  tj]  . 

dp 
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Note  that  since  E./E~  >>  1,  t(p)  -  t  must  be  negative  for  p  <  p  in 
1  *  1  0 

order  that  equation  (9)  reflect  a  frozen  A  B  reaction.  Impose  then 


t(P)  <  t1  for  p  <  pQ  , 

~  ~  2 

and  define  the  0{1)  constant  D_  as  D„  =  6  Dl  (hence  n-  =  2).  The 

2  2  2  2  2 

B  -*■  C  structure  equations  then  become 


(10) 


2  2 

p  <  pQ  :  =  0,  =  D2M~2z  exp(-y-Q2z-t2) ,  and 

dp  dp 


d  z  ~  -2 

p  >  pg  :  Y  =  0,  — “  =  DM  z  exp(-Q  z-t  ) 

dp2  " 


Matching  with  the  outer  solution  provides  the  conditions 


dy  dz 

TT  +  -Y  ,  —  ♦  -Z  .  and 

dp  — 00  dp  — 


dz 

p++o"sXT+0'  z  *  0 

dp 


Also  y,  z,  and  d(y  +  z)/dp  must  be  continuous  at  p^. 


y(P)  =  (P0  -  P)Y_ 


for  p  <  p 


0  ' 


z(pQ)  =  z(pQ)  =  z(PQ)  >  and 


dz  .  dz  .  +. 

dp  (PV  =  dp  (p0*  +  Y- 


It  follows  that 


The  value  of  p^  is  not  yet  known,  however.  It  must  be  fixed  at  the 

point  where  the  A  B  reaction  becomes  important,  i.e.  where 

T  -  T^  =  O(E^).  Thus,  t(p^)  ”  must  vanish.  This  fact  guarantees  (10) 
2  2 

since  d  t/dp  <  0  and  dt/dp  is  seen  to  be  positive  at  p^.  It  also  gives 
Q2z(P0>  =  Since  z  is  the  leading-order  term  in  Z,  z(pq)  must  be 

positive,  requiring  t^  to  be  negative,  which  is  not  surprising.  A 
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negative  requires  that  Tj  remain  less  than  through  0< ^ 2  > ' 

ensuring  that  the  A  ♦  B  reaction  remain  frozen  throughout  the  region 
preceding  its  zone,  as  is  necessary. 

The  B  ♦  C  structure  problem  may  now  be  written: 

for  P  <  PQ  "I 


i_z 

dp2 


r  D2M  2z  exp[  (p  -  Pq)*^  -  Q2z  -  t2] 


~  -2 

(  D2M  z  exp( -Q2 z  -  t2)  for  P  >  PQ 


z  +  z  p 

_ao 


"Y-po 


as  p 


z(p  )  —  — o  ^t 
u2  1 


t<v- 


dz 

dp 


<%> 


+  Y 


dz 

dp 


0,  z+0  as  p  ♦  +•  . 


(11) 


**»/  _ 

If  D2exp(-t2)  is  labeled  D2  and  t^  is  labeled  -t0,  (11)  takes  on  the 

form  of  the  B  +  C  structure  problem  obtained  and  numerically  solved  in 
[ 1 ] .  Those  numerical  results  give  the  equivalent  of 


M2  =  D2exp(-t2  )/F(-t1 ) 

=  E22T^D2exp(-E2/To)/F(-t1) 


(12) 


where  F  exhibits  the  asymptotic  behavior 


F  ~ 


4  Q2<Y-  +  Z-)2  “  *1  +  *“ 


as  t1  ♦  0 


(13) 


and  is  illustrated  in  the  graph  in  Figure  3  for  P2  =  .5  and  Y^  -  .75 
(Ref.  {1]  includes  an  extra  factor  of  ■j-  in  the  behavior  of  F  as  t^  ♦  0). 
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T  +  0(6.).)  It  is  therefore  possible  that  the  corresponding  burning 


rate  (8)  may  be  recovered  from  (12)  and  (13)  if  the  proper  distinguished 
limits  describing  t1  ♦  0  and  t2  *  -*•  are  considered.  This  requires  a 
careful  asymptotic  analysis,  however,  and  is  not  pursued  here.  It  may  be 
noted,  though,  that  the  limit  t^  0  gives  z(pg)  ♦  0  in  (11),  implying 
completion  of  the  B  ♦  C  reaction,  that  is  pg  ♦  “.  This  indicates  that 
the  A  ♦  B  flame  moves  to  the  right  of  the  B  C  flame. 

It  now  appears  feasible  that  either  of  the  two  possible  separated-flame 
configurations  may  be  recovered  from  the  merged-flame  configuration  just 
considered.  It  follows  that  it  may  be  possible  for  separated  flames  to  merge 
and  subsequently  separate  with  their  ordering  reversed;  i.e.  the  flames  may 
cross  as  T^  and  T2  are  varied.  This  conclusion  contradicts  that  of  Kapila 
and  Ludford. 

Although  the  result  (12)  has  the  identical  form  of  the  burning  rate  found 
in  [1],  there  is  one  crucial  difference  between  the  two  results:  the 
nonnegative  argument  of  F .  In  the  present  work,  the  argument  of  F  is 

"*1  =  *21(T»  '  V 


and  is  naturally  nonnegative.  The  argument  of  F  in  [1],  on  the  other  hand, 
is 


*0  ’  52,<T2  -  V  ■ 


so  that  the  restriction  tQ  >  0  leads  to  the  conclusion  that  T^  <  T2  when 
Ej  >>  E2,  admitting  only  one  type  of  separated  flames.  Indeed,  if  T1 
cannot  exceed  T2,  the  separated  flame  configuration  represented  by  (8) 
cannot  exist,  hence  the  flames  cannot  cross.  It  follows  that  the  ordering  of 
separated  flames  is  fixed  by  the  ordering  >>  E2 •  The  source  of  this 
restrictive  conclusion  is  Kapila  and  Ludford's  unspoken  assumption  that  T2 
exactly  equals  T—,  i.e.  that  t2  =  0. 
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The  structure  analysis  for  the  thinner  flame  contributes  only  information 
about  higher -order  terms  in  the  outer  solution.  Treatment  of  the  inner 
structure  problem  is  therefore  omitted. 

Analogous  results  and  conclusions  are  obtained  from  the  analysis  of 
merged  flames  with  E2  >>  E1 .  In  [1],  it  is  found  that  Tj  cannot  exceed 
T.|  in  that  case,  hence  only  separated  flames  represented  by  (8)  are 
possible.  This  follows  from  the  assumption  that  T^  is  exactly  T^.  If  a 
more  general  T^  is  allowed,  i.e.,  T1“Ta,  +  ^‘r^*  Bay,  the  restriction  on 
admissible  separated  flames  need  no  longer  be  true. 

Conclusion 

The  present  investigation,  of  which  the  analysis  in  [1]  constitutes  a 
special  case,  assumes  general  forms  for  the  characteristic  temperatures  T1 
and  T2  by  distinguishing  T^  (i  -  1,2)  from  its  limit  at  infinite 
activation  energies.  The  burning-rate  results  of  Kapila  and  Ludford  are  found 
to  be  true  in  general  for  separated  flames,  but  must  be  modified  for  merged 
flames.  Consequently,  some  of  their  conclusions  are  not  generally 
applicable.  In  particular,  their  analysis  of  merged  flames  forces  them  to 
conclude  that  the  order  in  which  separated  flames  occur  is  dictated  by  the 
relative  sizes  of  the  activation  energies.  This  restriction  is  a  consequence 
of  their  assumption  that  one  of  the  characteristic  temperatures  is  strictly 
equal  to  its  limit.  By  accounting  for  the  difference,  this  work  removes  that 
restriction. 
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